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Quantum dots may display fascinating features of strong correlation such as finite-size 
Wigner crystallization. We here review a few electron spectroscopies and predict that 
both inelastic light scattering and tunneling imaging experiments are able to capture 
clear signatures of crystallization. 
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1. Quantum Dots as Tunable Correlated Systems 

Semiconductor quantum dots^^^ (QDs) are nanostructures where electrons (or holes) 
are confined by electrostatic fields. The confinement field may be provided by com- 
positional design (e.g., the QD is formed by a small gap material embedded in a 
larger-gap matrix) or by gating an underlying two-dimensional (2D) electron gas. 
Different techniques lead to nanometer size QDs with different shapes and strengths 
of the confinement. As the typical de Broglie wavelength in semiconductors is of the 
order of 10 nm, nanometer confinement leads to a discrete energy spectrum, with 
energy splittings ranging from fractions to several tens of meV. 

The similarity between semiconductor QDs and natural atoms, ensuing from the 
discreteness of the energy spectrum, is often pointed out.'^-'*^^'^'^ Shell structure, 
correlation effects,^ Kondo physics, ^"^^^^ are among the most striking experimen- 
tal demonstrations. An intriguing feature of these artificial atoms is the possibil- 
ity of a fine control of a variety of parameters in the laboratory. The nature of 
ground and excited few-electron states has been shown to vary with artificially 
tunable quantities such as confinement potential, density, magnetic field, inter-dot 
coupling. ^'^■'''-'^^'^^'^^ Such flexibility allows for envisioning a vast range of applica- 
tions in optoelectronics (single-electron transistors,^^ lasers, micro-heaters and 
micro-refrigerators based on thermoelectric effects, ^^) life sciences, as well as 
in several quantum information processing schemes in the solid-state environment 
(e.g. Ref. 19). 
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Almost all QD-based applications rely on (or are influenced by) electronic corre- 
lation effects, which are prominent in these systems. The dominance of interaction 

in artificial atoms is evident from the multitude of strongly correlated few-electron 
states measured or predicted under different regimes: Fermi liquid, Wigner molecule 
(the precursor of Wigner crystal in 2D bulk), charge and spin density wave, incom- 
pressible state reminescent of fractional quantum Hall effect in 2D (for reviews see 
Refs. 12, 20). The origin of strong correlation effects in QDs is the following: While 
the kinetic energy term of the Hamiltonian scales as r~^, being the parameter 
measuring the average distance between electrons, the Coulomb energy term scales 
as rj^. Contrary to natural atoms, in QDs the ratio of Coulomb to kinetic energy 
can be rather large (even larger than one order of magnitude), the smaller the car- 
rier density the larger the ratio. This causes Coulomb correlation to severely mix 
many different Slater determinants. 

The way of tuning the strength of correlation in QDs we focus on in this paper is 
to dilute electron density. At low enough densities, electrons evolve from a "liquid" 
phase, where low-energy motion is equally controlled by kinetic and Coulomb energy, 
to a "crystallized" phase, reminescent of the Wigner crystal in the bulk, where 
electrons are localized in space and arrange themselves, in absence of disorder, in 
a geometrically ordered configuration (Wigner molecule^^), so that electrostatic 
repulsion is minimized^' So far, there are no direct experimental confirmations of 
the existence of these fascinating states. 

In this paper we focus on theoretical properties of electron states in Wigner 
molecules and on their possible signatures in experimentally available spectro- 
scopies. We predict that both inelastic light scattering and wave function imaging 
techniques, like scanning tunneling spectroscopy (STS), may provide direct access 
to the peculiar behavior of crystallized electrons. Specifically, the former is able 
to probe the "normal modes" of the molecules, while the latter is sensitive to the 
spatial order of the electron phase. 

The structure of the paper is as follows. After an exposition of our theoretical 
and computational approach (Sec. 2), we briefly review inelastic light scattering and 
imaging spectroscopies (Sec. 3), showing a few results for the four-electron Wigner 
molecule (Sees. 4 and 5). 

2. Theoretical Approach to the Interacting Problem: 
Configuration Interaction 

The theoretical understanding of QD electronic states in a vaste class of devices is 
based on the envelope function and effective mass model. Here, changes in the 
Bloch states, the eigenhmctions of the bulk semiconductor, brought about by "ex- 
ternal" potentials other than the perfect crystal potential, are taken into account 
by a slowly varying (envelope) function which multiplies the fast oscillating peri- 
odic part of the Bloch states. This decoupling of fast and slow Fourier components 
of the wave function is valid provided the modulation of the external potential is 
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slow on the scale of the lattice constant. Then, the theory allows for calculating 
such envelope functions from an effective Schrodinger equation where only the ex- 
ternal potentials appear, while the unperturbed crystalline potential enters as a 
renormalized electron mass, i.e. the effective value m* replaces the free electron 
mass. Therefore, single-particle (SP) states can be calculated in a straightforward 
way once compositional and geometrical parameters are known. This approach was 
proved to be remarkably accurate by spectroscopy experiments for weakly confined 
QDs;^^'^^ for strongly confined systems, such as certain classes of self-assembled 
QDs, atomistic methods might be necessary.^^ 

The QD "external" confinement potential originates either from band mismatch 
or from the self-consistent field due to doping charges. In both cases, the total field 
modulation which confines a few electrons is smooth and can often be approximated 
by a parabolic potential in two dimensions. Then, the fully interacting effective-mass 
Hamiltonian of the QD system reads as:^^ 

with 

^oW = ^ + ^™*^oV?. (2) 

Here, A'^ is the number of free conduction band electrons localized in the QD, e and 
K are respectively the electron charge and static relative dielectric constant of the 
host semiconductor, r is the position of the electron, p is its canonically conjugated 
momentum, loq is the natural frequency of a 2D harmonic trap. 

The eigenstates of the SP Hamiltonian (2) are known as Fock-Darwin (FD) 
orbitals.^ The typical QD lateral extension is given by the characteristic dot radius 
^QD = {h/m*uJoY/'^, Iqt) being the mean square root of r on the FD lowest-energy 
level. As we keep A'^ fixed and increase ^qd (decrease the density), the Coulomb-to- 
kinetic energy dimensionless ratio^^ A = ^qd/cib [O'b = ti?K,/{m*e^) is the effective 
Bohr radius of the dot] increases as well, driving the system into the Wigner regime. 

Wo solve numerically the few-body problem of Eq. (1), for the ground (or ex- 
cited) state at different numbers of electrons, by means of the configuration interac- 
tion (CI) method:^'* We expand the many-body wave function in a series of Slater 
determinants built by filling in the FD orbitals with N electrons, and consistently 
with global symmetry constraints. Specifically, the global quantum numbers of the 
few-body system are the total angular momentum in the direction perpendicular 
to the plane, M, the total spin, 5, and its projection along the 2;-axis, Sz- In the 
Slater-determinant basis the few-body Hamiltonian (1) is a large;, block diagonal 
sparse matrix that we diagonalizc by means of a newly developed parallel code.^'"' 
Note that, before diagonalization, we are able to build separate sectors of the Fock 
space corresponding to different values of {M,S,Sz). We use, as a single-particle 
basis, up to 36 FD orbitals, and we are able to diagonalize matrices of linear di- 
mensions up to « 10^ (see Ref. 24). The output of calculations consists of both 
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energies and wave functions of the selected correlated states. We then post-process 
wave functions to obtain those response functions connected to the spectroscopy of 
interest, which we will discuss in Sec. 3. 

3. Quantum Dot Spectroscopies 

There are two main classes of electron spectroscopies in QDs (cf. Fig. 1). In single- 
electron tunneling spectroscopies one is able to inject just one electron into the 
interacting system, in virtue of the Coulomb blockade effect. In this way one 
accesses quantities related mainly to the ground states of both N and -|- 1 elec- 
trons, like the dot chemical potential Eq{N -|- 1) — Eq{N), where Eo{N) is the 
A''-electron ground state energy. While such transport spectroscopies were very suc- 



transport Eo(N+1) - Eo(N) 




optics E'(N) - Eo(N} 




Fig. 1. Transport vs. optical electron spectroscopies in quantum dots. 

cessful in QD studies, they suffer severe limitations when probing excited states. On 
the other hand, optical techniques such as far-infrared and inelastic light scattering 
spectroscopies, which leave the probed system uncharged, allow for an easier access 
to excitation energies E*{N). Below we focus on two specific examples of these 
families of experiments (for reviews see Refs. 1, 2, 12, 20). 

3.1. Inelastic light scattering 

Inelastic light scattering experiments in QDs probe low-lying neutral excitations. 
Depending on the relative orientation of the polarizations of the incoming and 
scattered photons, one is able to access either charge or spin density fluctuations.^^ 
Different excited states can be probed by varying the frequency of the laser with 
respect to the optical energy band gap of the host semiconductor, or by changing 
the momentum tranferred from photons to electrons. ^'^^'^^'^^ 

Only recently QDs with very few electrons were studied. ^'^^ Here we focus on 
charge density excitations probed when the laser frequency is far from the resonance 
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with the optical gap.^°'^^ In this hmit, the scattering cross section da/duj at a given 
energy w is proportional to the momentum-resolved dynamical response function: 

^ « ^ |M„o|'(5(a; +a;o), with M„o = Zaba ("le'^'l^) (n|ct^Cb,|0) . (3) 

n 

Here |0) and |n) are the ground and excited interacting few-body states, as obtained 
by the CI computation, with energies wq and w„, respectively, cl^ is a fermionic 
operator creating an electron in the o-th FD orbital with spin a, q is the wave 
vector transferred in the inelastic photon scattering event. Formula (3) describes 
how ground and excited states are selectively coupled by charge density fluctuations 
of momentum q. 



3.2. Wave function imaging via tunneling spectroscopy 

The imaging experiments, in their essence, measure quantities directly proportional 
to the probability for transfer of an electron through a barrier, from an emitter, 
where electrons fill in a Fermi sea, to a dot, with completely discrete energy spec- 
trum. In multi-terminal setups one can neglect the role of electrodes other than the 
emitter, to a first approximation. The measured quantity can be the current, ^^'^^ 
the differential conductance, ^'^^'^^ or the QD capacitance, ^'^^'^^ while the emitter 
can be the STS tip,^^-'''*-''''^ or a n-doped GaAs contact,"^-'^'-'^-'^'^'^'^^ and the barrier 
can be the vacuum"^^''^^''^^ as well as a AlGaAs spacer. ^'^'•^■^'•^^''^'^ 

The measured quantities are generally understood as proportional to the den- 
sity of carrier states at the resonant tunneling (Fermi) energy, resolved in either 
reaP^''^^'*^^ or reciprocal"^^'"^^ space. However, Coulomb blockade phenomena and 
strong inter-carrier correlation — the fingerprints of QD physics — complicate the 
above simple picture. Below we review the theoretical framework we have recently 
developed"'^ '■'^ to clarify the quantity actually probed by imaging spectroscopies. 

According to the seminal paper by Bardeen,*^ the transition probability (at 
zero temperature) is given by the expression (27r//i) \M\ n{ef), where M. is the 
matrix element and n(e/) is the energy density of the final QD states. The standard 
theory would predict the probability to be proportional to the total density of QD 
states at the resonant tunneling energy, e/, possibly space-resolved since M would 
depend on the resonant QD orbital.^^ Let us now assume that: (i) Electrons in the 
emitter do not interact and their energy levels form a continuum, (ii) Electrons 
from the emitter access through the barrier a single QD at a sharp resonant energy, 
corresponding to a well defined interacting QD state, (iii) The QD is quasi-2D, the 
electron motion being separable in the xy plane and z axis, which is parallel to the 
tunneling direction, (iv) Electrons in the QD all occupy the same confined orbital 
along z, Xqb{z)- Then one can show^^ that the matrix element Ad may be factorized 
as 



M<xTM, 



(4) 
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where T is a purely single-particle matrix element while the integral M contains 
the whole correlation physics. 

The former term is proportional to the current density evaluated at any point 
2bar in the barrier: 

(5) 

Z = 2bar 

where Xe(^) is the resonating emitter state along z evanescent in the barrier. The 
term (5) contains the information regarding the overlap between emitter and QD 

orbital tails in the barrier, xe{z) and xqd{z), respectively. Since T is substantially 
independent from both TV and xy location, its value is irrelevant in the present 
context. 

On the other hand, the in-plane matrix element M conveys the information 
related to correlation effects. If we now specialize to STS and assume an ideal, 
point-like tip, then M oc (fiQuif) where </'Qd(t") is the quasi-particle (QP) wave 
function of the interacting QD system: 

<^QD(r) = (A'--l|*(r)|iV). (6) 

Here \I/(r) is the fermionic field operator destroying an electron at position r = 
{x,y), I A'' — 1) and |A'') are the QD interacting ground states with N — 1 and N 
electrons, respectively, calculated via CI (Sec. 2). We omit spin indices for the sake 
of simplicity. 

Therefore, the STS differential conductance is proportional to \ipQ-£,{r)\'^ , which 
is the usual result of the standard non-interacting (or mean-field) thcory'^''''"'^ for 
the highest-energy occupied (or lowest-energy unoccupied) orbital, provided the SP 
(mean-field) orbital is replaced by the QP wave function unambiguously defined by 
Eq. (6). 

4. Inelastic Light Scattering Spectroscopy in the Wigner Limit 

We focus on the four-electron system at density low enough (A = 10) to induce 
the crystallization of a Wigner molecule. Figure 2(a) displays our CI results for the 
low-energy region of the excitation spectrum for different values of the total orbital 
angular momentum and spin multiplicities (the so called yrast line^^). The ground 
state is the spin triplet state with M = 0, which is found to be the lowest-energy 
state in the range < A < 20. The absence of any level crossing as the density 
is progressively diluted (i.e. A increases) implies that the crystallization process 
evolves in a continuos manner, consistently with the finite-size character of the 
system. Nevertheless, several features of Fig. 2(a) demonstrate the formation of a 
Wigner molecule in the dot. 

First, the two possible spin multiplets other than 5=1, namely S = 0,2, lie 
very close in energy to the ground state. This is understood by invoking electron 
crystallization. In fact, in the Wigner limit, the Hamiltonian of the system turns 



2m* 



dz 



dz 
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Fig. 2. (a) Excitation energies of tlic four-electron Wigncr molecule vs. total orbital angular 
momentum M, for different values of the total spin S, with A = 10. (b) Corresponding Raman 
scattering cross section for charge excitations in the ofT-resonance limit. Here we set |q| = 2 x 10* 



cm ^ and introduce a fictitious gaussian broadening of Raman peaJcs (u : 
parameters were used throughout the paper. 



: 0.9 X IQ-^ eV). GaAs 



into a classical quantity, since the kinetic energy term in it may be neglected with 
respect to the Coulomb term. Therefore, only commuting operators (the electron 
positions) appear in the Hamiltonian. In this regime the spin, which has no classical 
counterpart, becomes irrelevant:^^'^^ spin-dependent energies show a tendency to 
degeneracy. This can also be understood in the following way: if electrons sit at 
some lattice sites with unsubstatial overlap of their localized wave functions, then 
the total energy must not depend on the relative orientation of neighboring spins. 

Second, an exact replica of the sequence of the two lowest-energy singlet and 
triplet states occurs for both M = and for M = 4 [(cf. the triplet state labeled A 
in Fig. 2(a)]. Such a period of four units on the M axis identifies a magic number, ^'^^ 
whose origin is brought about by the internal spatial symmetry of the interacting 
wave function. ^^'^^ In fact, when electrons form a stable Wigncr molecule, they 
arrange themselves into a four-fold symmetry configuration, where charges are lo- 
calized at the corners of a square.^^'*^'*^ 

A third distinctive signature of crystallization is the appearance of a rotational 
band,'^^ which in Fig. 2(a) is composed of those lowest-energy levels, separated 
by energy gaps of about 0.1 meV, that increase monotonically as M increases. 
This band is called "rotational" since it can be identified with the quantized levels 
Erot{M) of a rigid two-dimensional top, given by the formula 
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where / is the moment of inertia of the top.*** These excitations may be thought of 
as the "normal modes" of the Wigner molecule rotating as a whole in the xy plane 
around the vertical symmetry axis parallel to z. 

Figure 2(b) displays the calculated off-resonance Raman spectrum for charge 
density fluctuations {AS = 0, where AS is the variation of the spin with respect 
to the ground state). The dominant peak is the one labeled A, corresponding to 
the AM — 4 normal mode of the Wigner molecule of Fig. 2(a), which therefore 
represents a clear, visible feature of Wigner crystallization. Notice also the appear- 
ance of the so called Kohn mode^ (labeled K) at higher energy, which is a dipolar 
(AM = 1) collective motion of the center of mass of the electron system. 

5. Wigner Molecule Formation Seen by Imaging Spectroscopy 

Figure 3 shows the QP wave function, corresponding to the injection of the fourth 
electron into the QD, as a function of x {at y = 0), for two different values of A. As 
A increases (from 0.5 to 10), the density decreases going from the non- interacting 
limit (A = 0.5) deep into the Wigner regime (A = 10). At high density (A = 0.5, 
approximately corresponding to the electron density Ug = 3.3 x 10^^ cm~^) the wave 
function substantially coincides with the non-interacting FD p-like orbital. By in- 




-4-2 2 



^ Cqd) 

Fig. 3. Quasi-paiticlc wave function vs. x (y = 0) for different values of the dirnensionless 
Coulomb-to-kinetic energy ratio A. The STS differential conductance, proportional to the wave 
function square modulus, corresponds to the tunneling process N = 3 —* N = 4. The length unit 
is the characteristic dot radius ^qd. 

creasing the QD radius (and A), the QP wave function weight moves towards larger 
values of r (the s-like symmetry of the A = 10 curve comes from a transition of the 
three-electron ground statc^'* around A « 6). By measuring lengths in units of £qd, 
as it is done in Fig. 3, this trivial effect should be totally compensated. However, 
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we see for A = 10 (rig ~ 1.1 x 10^ cm^^) that the now much stronger correlation 
is responsible for an unexpected weight reorganization, which is related to the for- 
mation of a "ring" of crystallized electrons in the Wigncr molecule (cf. Sec. 4). The 
shape of the A = 10 curve of Fig. 3 is consistent with the onset of a solid phase 
with four electrons sitting at the apices of a square, as discussed in Sec. 4. Note also 
the dramatic weight loss of QPs as A is increased: the stronger the correlation, the 
more effective the "orthogonality" between interacting states. 

6. Conclusions 

The Wigner molecule is an intriguing electron phase, peculiar of QDs, which still 
lacks for experimental confirmation. We predict that both inelastic light scattering 
and imaging spectroscopies are able to probe distinctive features of crystallization. 
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